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Abstract 

We investigate the role weakly damped modes play in the selection of Faraday wave 
patterns forced with rationally-related frequency components rmo and noj. We use 
symmetry considerations to argue for the special importance of the weakly damped 
modes oscillating with twice the frequency of the critical mode, and those oscillat- 
ing primarily with the "difference frequency" \n — m|u; and the "sum frequency" 
(n+m)u). We then perform a weakly nonlinear analysis using equations of Zhang and 
Vihals [1] which apply to small-amplitude waves on weakly inviscid, semi-infinite 
fluid layers. For weak damping and forcing and one-dimensional waves, we perform 
a perturbation expansion through fourth order which yields analytical expressions 
for onset parameters and the cubic bifurcation coefficient that determines wave am- 
plitude as a function of forcing near onset. For stronger damping and forcing we 
numerically compute these same parameters, as well as the cubic cross-coupling 
coefficient for competing waves travelling at an angle 9 relative to each other. The 
resonance effects predicted by symmetry are borne out in the perturbation results 
for one spatial dimension, and are supported by the numerical results in two dimen- 
sions. The difference frequency resonance plays a key role in stabilizing superlattice 
patterns of the SL-I type observed by Kudrolli, Pier and Gollub [2] . 
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1 Introduction 



The Faraday wave system provides the canonical example of how spatiotem- 
poral patterns form through a parametric instability. In this system, a fluid 
subjected to a time-periodic vertical acceleration of sufficient strength under- 
goes an instability to standing waves on the free surface. In his original exper- 
iment [3] Faraday observed that the standing waves had half the frequency of 
the forcing; this is the familiar subharmonic response. Other experimentalists 
subsequently observed familiar patterns such as stripes, squares, and hexagons 
(see [4] for a review). 

More recent experiments have utilized the two-frequency forcing function, 
which we may write in the following forms: 

g{t) = gz[cos{x) cos^mut) + sm{x){n^t + 4>)] (1) 




Here m and n are co-prime integers, so that the forcing function is periodic 
with period T = 2ti/uj. An interesting feature of this forcing function is that 
the primary instability leading to Faraday waves may be either harmonic or 
subharmonic (with respect to T) depending on the value of % and the parities 
of m and n. For instance, if the cos(ma;t) component is dominant and if m 
is even (odd) then the bifurcation will be to harmonic (subharmonic) waves. 
This was demonstrated numerically by the linear stability analysis of Besson 
et al. [5]. For m and n not both odd, there is a codimcnsion-two point in the 
Qz - X parameter space (or alternatively, in Qm - Qn space) at which harmonic 
and subharmonic instabilities occur simultaneously at different spatial wave 
numbers. The corresponding value x = Xbc is called the "bicritical point" . Ex- 
periments performed near the bicritical point have produced exotic patterns, 
including triangles [6], quasipatterns [2,7] and superlattice patterns [2,8-10]. 

The term "superlattice pattern" refers to a periodic pattern that has spatial 
structure on more than one length scale - typically, a small scale structure 
and a large scale spatial periodicity. A wealth of recent experimental work 
has produced superlattice patterns. In addition to the Faraday experiments 
with two-frequency forcing mentioned above, superlattice patterns have been 
observed in nonlinear optical systems [11], in vertically vibrated Rayleigh- 
Benard convection [12,13], in Faraday waves with single frequency forcing at 
very low frequencies [14] and in granular layers forced with two frequencies [15]. 
The superlattice patterns come in a variety of types, and may be comprised of 
different numbers of critical modes having different spatial arrangement and 
temporal dependence. 




Qm cos{mujt) + gn{nujt + 0) 
G^e^'"^* +G'„e*"'^*+c.c. 
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We mention two types of superlattice patterns here. One type is called "su- 
perlattice two" or SL-II in [2]. Patterns of this type have been shown to arise 
in a spatial period-multiplying bifurcation from hexagons; see [16,17] for a 
bifurcation analysis of these patterns. In contrast, "superlattice one" or SL-I 
patterns exist as primary branches bifurcating directly from the trivial so- 
lution [18]. For SL-I Faraday wave patterns, one of the length scales is set 
(approximately) by the dominant frequency component in (1). Curiously, the 
second length scale in the pattern is not set by the other forcing frequency. 
For instance, a numerical linear analysis using the parameters corresponding 
to the experimental SL-I pattern from [2] indicates that the two critical wave 
numbers are in a ratio of 1.22 at the bicritical point, but the ratio of the two 
length scales in the pattern is actually observed to be a/7 ~ 2.65. The iden- 
tification of a mechanism for the selection of the second length scale is our 
motivation for this paper. 

Many studies of Faraday wave pattern formation [1,6-10,14,19-23] focus on 
the effect that resonant triad interactions have on the nonlinear pattern se- 
lection problem. The triads of spatially resonant modes satisfy ki ± k2 = ks, 
where |ki| = |k2| is the wave number of one of the instabilities at the bicritical 
point, and Iks] is the wave number of the other. In [24] temporal symmetry 
arguments were used to show that this triad interaction affects pattern selec- 
tion on the subharmonic side of the bicritical point, but not on the harmonic 
side. In [25] it was shown that weakly damped harmonic modes not associ- 
ated with the bicritical point can affect harmonic and subharmonic pattern 
selection. Explicit numerical calculations were performed to demonstrate that 
the presence of a weakly damped harmonic mode influences pattern selection 
by affecting the cubic cross-coupling coefficient in the bifurcation equations 
describing the dynamics of competing waves. 

In this paper, too, our goal is to investigate the role weakly damped modes 
play in pattern selection. We use symmetry considerations to explain the spe- 
cial importance of particular weakly damped harmonic modes in terms of 
their contributions to cubic cross-coupling coefficients in the relevant bifurca- 
tion equations. Specifically, we find that the most important weakly damped 
modes are the temporal harmonic of the Faraday-unstable mode oscillating 
with dominant frequency muj, the "difference frequency mode" oscillating with 
dominant frequency \m — n\uj^ and the "sum frequency mode" oscillating with 
dominant frequency {m+n)uj. (Here, we have assumed without loss of general- 
ity that the cos{mut) component in (1) is the dominant one. Additionally, we 
exclude the forcing frequency ratios m/n or n/m equal to 1/2, 1/3, 1/4, and 
2/3, in which case strong resonances are possible). Our symmetry arguments 
are borne out in a weakly nonlinear analysis that uses equations of Zhang and 
Vihals [1] which apply to weakly inviscid, semi-infinite fiuid layers. For weak 
damping and forcing and one- dimensional waves, we perform a perturbation 
expansion through fourth order which yields analytical expressions for onset 
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parameters and the cubic self-interaction coefficient that determines wave am- 
phtude as a function of forcing amphtude near onset. For stronger damping 
and forcing we compute these same parameters numerically as well as the cu- 
bic cross-coupling coefficient B{9) for competing waves travelling at an angle 
9 relative to each other. Prom the resulting analytical expressions for the one- 
dimensional case, we are able to quantify the effect of the key resonances and 
see how their existence depends on the forcing frequency ratio m/n. For the 
two dimensional case, our numerical results show that the resonance effects 
follow the same scaling laws as in the 1-d case. A simple argument, vahd for 
weak damping and forcing which relies only on the inviscid dispersion rela- 
tion, allows us to predict the spatial angles at which the resonances occur, 
and to see how their existence depends on m, n and a dimensionless fluid 
gravity-capillarity parameter. A bifurcation analysis reveals that the differ- 
ence frequency resonance can help stabilize an SL-I pattern whose large-scale 
periodicity depends on the wavelength associated with the difference frequency 
mode. We note that the difference frequency mode has been observed to play 
an important role in experiments where other complex Faraday wave patterns 
are formed [8-10]. 

This paper is organized as follows. In section 2.1 we review basic ideas about 
resonant triad interactions and their potential for affecting SL-I pattern se- 
lection. In section 2.2 we use symmetry arguments to identify which weakly 
damped modes we expect to be most important in terms of their contributions 
to the cross-coupling coefficient B{9). Section 3 contains the weakly nonlin- 
ear analysis of the Zhang- Vinals equations for weak damping and forcing and 
one-dimensional waves, which leads to approximate formulas for the critical 
forcing and wave number, and the cubic self-interaction coefficient. The per- 
turbation results for onset parameters are discussed in 4.1. The expression 
for the self-interaction coefficient and numerical results for the cross-coupling 
coefficient are examined in sections 4.2.1 and 4.2.2 respectively, with special 
attention given to the role played by resonant triads and the implications for 
SL-I pattern selection. We summarize our main results in section 5. 



2 Background 

2.1 Resonant triads, standing wave equations and pattern stability 

For Faraday waves on a domain of infinite horizontal extent there is no pre- 
ferred direction, so each wave number is associated with a circle of wave vectors 
in Fourier space. One of the simplest mechanisms through which waves on two 
different Fourier circles may interact is a resonant triad interaction. Such res- 
onant triads consist of three wave vectors that determine an associated angle 
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Fig. 1. Spatially resonant triads of three wave vectors and their associated resonant 
angle Or- (a) ki + k2 = ks. (b) ki - k2 = ka. 

^.e(0, f]. 

Two examples of spatially resonant triads are shown in figure 1. In figure la, 
the spatial resonance condition is 

ki + k2 = ks. (2) 



Here, ki and k2 are wave vectors associated with neutrally stable modes of the 
Faraday instability, so that |ki| = |k2| = /cc- The wave vector [ksl is associated 
with a damped harmonic mode. The resonant angle Or satisfies 

(3) 




2|ki| 



For this resonant triad, v/2|ki| < jkal < 2. It is also possible to have a resonant 
triad where < jkal < \/2|ki|; see figure lb. Here, the spatial resonance 
condition is 

ki - k2 = kg. (4) 



and the resonant angle 9r satisfies 




We follow [24,25] and consider equations describing the slowly-varying ampli- 
tudes of the modes with wave vectors ki . . . ka, which we assume satisfy the 
spatial resonance condition (2) (the case for (4) is similar): 
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= AlZi + q;iZ2^3 + {A\Zi\'' + + C\Z3\^)Zi (6) 

= A1Z2 + aiZiZa + {A\Z2\^ + h\Z^f + (71^31^)^2 
Z3 = A2Z3 + aiZ^Z-i + (L>|Zi|2 + D\Z2f + £;|Z3|^)Z3. 

Here Zi, Z2, and ^3 are the slowly varying amplitudes of the Fourier modes 
with wave vectors ki, k2, and k3 respectively, and the dot represents differ- 
entiation with respect to the slow time scale. All coefficients are real-valued. 
If the Z3 mode is in fact damped (i.e. A2 < 0) and the Zi and Z2 modes are 
neutrally stable (i.e. Ai = 0) then a further center manifold reduction may be 
performed to the critical Z\ and Z2 modes. Then Z^ satisfies 

^3 = -^^1^2 + ..., (7) 
A2 

and the (unfolded) bifurcation problem, to cubic order, is 

" =AiZi + A|Zi|2Zi + S(^,)|Z2|2Zi (8) 



dT 
dZ2 

It 

where 



AlZ2 + A\Z2\^Z2 + B{9r)\Zi\^Z2, 



B{9r) = h-^. (9) 



The dependence on the resonant angle indicates that the cross-coupling 
coefficient B{9) for competing waves travelling at an angle 9 relative to each 
other is evaluated at the angle of spatial resonance. 

For appropriately chosen angles 9^., the bifurcation problem (8) describes dy- 
namics on an invariant subspace of the twelve-dimensional problem which de- 
scribes the competition of simple rolls, simple hexagons, certain rhombic pat- 
terns, and SL-I patterns. We are ultimately concerned with the cross-couphng 
coefficient B{9), which plays a key role in determining the relative stabihty 
of these patterns. In particular, the stability properties of SL-I patterns with 
characteristic angles 9^ are enhanced when B{9h) is small in magnitude; see 
[25] for a detailed discussion. 

In [26], symmetry arguments are used to derive a scaling law for the quadratic 
coefficients ai and «2 in (6) when the resonant triad applies to the bicritical 
point, and for the case of weak damping and forcing. For instance, it is shown 
that for m odd and n even, and for the case that the cos(mcjt) forcing frequency 

n — 2 m — 1 

component dominates, ai and a2 are each proportional to gm Qn ^ . Thus 
depending on m and n, the quadratic terms in (6) can be quite small. Its 
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contribution to the pattern selection problem can only be made significant by 
getting sufficiently close to the bicritical point where A2 ^ 0. 

Here we will focus on resonant triads other than those associated with the 
bicritical point. In particular, we will identify resonant triads for which the 
quadratic terms scale (at most) linearly with g, i.e. the scaling is independent 
of m and n. Thus, we expect these triads to play a more important role in 
pattern formation, at least away from the bicritical point and for weakly forced 
waves. 



2.2 Determination of important resonances for weak damping and forcing 

Our goal in this section is to examine resonant triads from a symmetry per- 
spective with a special emphasis on temporal symmetries. Without loss of 
generality, we assume (unless otherwise specified) that the cos{mu!t) forcing is 
of greater significance than the cos{nujt) forcing. For weak damping, the crit- 
ical Faraday waves oscillate with a dominant frequency component of muj/2. 
We also consider weakly damped waves of frequency fl > {fl ^ muj/2), to 
be determined such that they lead to the largest possible contribution to the 
cross-coupling coefficient B{9) when slaved away; cf (7) - (9) in section 2.1. 
Strong resonances arc possible when the forcing frequency ratio m/n or n/m 
is equal to 1/2, 1/3, 1/4, and 2/3; we exclude these cases from our analysis. 

We follow [26] and focus on travelling waves, for which the action of time- 
translation is transparent. The travelling wave bifurcation equations are then 
reduced to those describing the standing wave problem. Specifically, we expand 
the fluid surface height /i(x, i), x e M^, in terms of the following six travelling 

waves: 

gi(ki-x+ima;t) gi(ki-x-ima;t) gi(k2-x+ima;t) ^-|^q^ 

+ W2 e*(^2-x-ima;t) gi(k3-x+nt) gi(k3-x-m) ^ ^ ^ 

Here, Zj and Wj, j = 1,2,3, are the slowly varying amplitudes of the travelling 
waves. The wave vectors ki . . . ks are assumed to satisfy the spatial resonance 
condition (2) (the argument for (4) is similar). The frequency Q and the wave 
number jkaj are related by a dispersion relation. In writing (10) we have as- 
sumed the problem is posed on an unbounded horizontal domain and then 
restricted our attention to solutions that are periodic on a rhombic lattice. 
Spatial translation symmetry acts on {zj, Wj), j = 1, 2, 3, as 

T(ei,e2) : {z„Wj) ^ {zj,wj)e'^^, 63 = 01 + 62 (11) 
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where (©i, ©2) £ T^. A rotation by tt, denoted by R, acts as 



TZ : {zj,Wj) {wj, Zj), j = 1, 2, 3. 



(12) 



and a reflection in the plane containing ka, denoted by k, acts as 



(13) 



(2:3, W3) {Z3,W3). 

Furthermore, there is a time-translation symmetry which acts on the forcing 
parameters in (1) and on the complex travelling wave amplitudes in (10): 



We now determine which quadratic terms will be allowed in the travelling wave 
amplitude equations, anticipating that these terms will lead to contributions 
to B{6). For example, from the spatial translation symmetry (11), the only 
quadratic terms that are allowed in the ii equation are Z2Z3, Z2W2,, W2Z3 and 

W2W3. 

We now consider the restrictions placed by the time translation symmetry (14) 
in order to determine which Q are allowed. We expect the largest contributions 
to B{6) in (8) to occur when the coefficients of quadratic terms are independent 
of the forcing amplitudes Gm and Gn at leading order, at least for small forcing. 
In this case, there is only one quadratic term that is permitted, namely 

i. W2W3 with Q — mw. 

The next largest contributions to B{9) occur when the coefficients of the 
quadratic terms in (6) are proportional to one power of G„i or Gn at leading 
order. In this case, the permitted equivariant terms in the Zi equation are 

ii. GmZ2Z3 with fl = 2muj 

iiia. GnZ2Z3 with fl = {m — n)u! if m > n 

iiib. GnZ2W3 with fl — {n — m)u) if n > m 
iv. GnZ2Zz with fl — {m-\- n)uj 

V. GnW2Z3 and GnW2W3 with D, — nuj. 

We may immediately dispense with several of these cases. The resonance in 
case ii. is not relevant for our investigation of Faraday waves because the 



Tm:{zuZ2,Z3)^{z^S'''' 

{Wi,W2,W3) (wie"*2 
{Gm, Gn) {Gm G**"' 




(14) 
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weakly damped mode oscillating with frequency 2muj is at sufficiently high 
wave number that the spatial resonance condition (2) cannot be satisfied for 
the inviscid dispersion relation. The resonance in case v. does not result in 
a contribution to B{9) at linear order in Gm, Gn- This may be understood 
by considering the effects of an approximate time reversal symmetry and an 
approximate Hamiltonian structure [27] and has been verified by an explicit 
perturbation calculation similar to those performed in section 3. 

We refer to iiia. and iiib. as cases of "difference frequency resonance." We now 
examine the amplitude equations for case iiia. (case iiib. is analogous) which 
are determined by the symmetries (11) - (14). The cubic truncation takes the 
form 



ii = XiZi + 5GmWi + Gn(^lZ2Zz + roWiZ2W2 (15) 

+ + r2\z2\^ + rsl^al^ + r^\wi\^ + rr,\w2\^ + rel^sD^i 

4 = A2^3 + GnP2ZlZ2 (16) 
+ (^^71^1^ + ^71^2^ + ^sl^sT + ^9^1^ + ^^91^2^ + ^loksH^S- 

Related equations for £2, Wi, W2, and W3 can be obtained from the discrete 
spatial symmetries (12) and (13). The SGmWi term in (15) is the usual para- 
metric forcing term. We have dropped linear and quadratic terms that scale 
higher than linearly in Gm and Gn- We have also dropped any cubic terms 
whose coefficients depend on these parameters. 

Since the resonant Z3 mode is damped (A2 < 0), we may slave it away so that 
(15) becomes 

zi = Ai^i + 5GmWi + roWiZ2W2 (17) 

|G„p/3l/32\ I |2 I I |2 I I 12I 

r2 ^ I F2I +r^\wi\ +r5|w2| > ^1. 

and equations related by the discrete spatial symmetries (12) and (13). 

For sufficiently large forcing \Gm\, the trivial solution of (17) loses stability. 
A center manifold reduction to standing waves equations of the form (8) may 
be formed at the critical forcing strength. The cross- coupling coefficient B{9) 
in (8) then includes a contribution proportional to |G„p/(ReA2) that results 
from slaving the difference frequency mode. 

Similar arguments can be made for case iv., in which the resonant mode oscil- 
lates at the so-called "sum frequency" (m + n)uj. For this case, too, the slaved 
mode results in a contribution to B{9) that is proportional to |G„p/(ReA), 
where Re A represents the damping of the slaved resonant mode. 
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Case i. corresponds to the well-known 1:2 temporal resonance, which is present 
for single frequency forcing [1]. An analysis similar to that performed above 
reveals that slaving of the damped mode oscillating with dominant frequency 
muj results in a contribution to B{9) which is independent of and Gn, and 
is inversely proportional to Re(A). Thus we expect that this contribution will 
be larger than that due to the sum or difference frequency resonances. 

For weak damping and forcing, then, we expect the standing wave modes 
with frequency mu, \m — n\uj, and (m + n)u! to be the most important weakly 
damped modes in terms of their contributions to B{9). Due to the constraints 
imposed by temporal symmetries, any other resonant modes will necessarily 
have higher powers of g in front of the quadratic terms in their travelling 
wave equations, and thus will result in smaller contributions to B{9) for weak 
forcing. The exception is the wave oscillating with dominant frequency |?^a;. 
Because this bifurcates directly from the bicritical point, its damping may be 
made arbitrarily small for appropriately chosen parameters, and its slaving 
can result in a large contribution to B{6) as demonstrated in [25]. Since our 
analysis assumes that the resonant modes have finite damping (i.e. we are 
bounded away from the bicritical point), this case is excluded here. 



3 Perturbation analysis for one dimensional waves 



Using the ideas discussed in the previous section, we now perform a per- 
turbation analysis on the Zhang- Vifials Faraday wave equations (which we 
introduce in section 3.1) to obtain quantitative results for Faraday waves in 
one spatial dimension. The forcing frequency ratios m/n ot n/m equal to 1/2, 
1/3, 1/4, and 2/3 are special cases of strong resonance which we exclude from 
our analysis. 

Since there is no spatial angle 6 to tune in one dimension, the 1:2 spatial 
resonance is the only possibility for the modes depicted in figure 1. For the 
one-dimensional case, then, a resonant triad interaction occurs when a stand- 
ing wave with critical wave number k and its spatial harmonic with wave 
number 2k fulfill one of the temporal resonance conditions from section 2.2. 
This situation may be achieved by varying fiuid parameters in the dispersion 
relation. In particular, in our calculations we vary a dimensionless capillar- 
ity parameter. Near those special values of the capillarity parameter where a 
temporal resonance occurs, we expect additional contributions to the cubic 
self-interaction coefficient A in the standing wave bifurcation equation 

^^A,Z, + A\Z^\'Z„ (18) 
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which is simply (8) restricted to one spatial dimension. 



3.1 The Zhang- Vinals Hydrodynamic Equations 



Zhang and Vinals derive from the Navier Stokes equations reduced equations 
for Faraday waves in [1]. This derivation is accomphshed by focusing on fluids 
of low viscosity and making a quasipotential approximation. The resulting 
equations apply to weakly damped, small amplitude surface waves on a semi- 
infinite layer of fluid. The system consists of two evolution equations for the 
surface height /i(x, r) and surface velocity potential $(x, r), where x e is 
the horizontal coordinate. We assume periodic boundary conditions. 

The Zhang- Viiials equations are 



(dr--fV^)h-V^^J^(h,^) (19) 
(dr - 7V^)$ - (FoV^ - G{t)) h = g{h, $), (20) 



where the nonlinear terms are 



:r(/i,$) = -V • (W$) + ^V2(/i^P$) -P(/iP$) (21) 
g{h, $) = -(D^f - -{V^y - (P$) {/iV^$ + V{hV^)} (22) 

-lroV-{{Vh){vhr}. 

(For brevity, we drop the h and $ dependence of and Q from now on.) The 
operator X> is a nonlocal operator that multiplies each Fourier component of 
a field by its wave number, e.g. T) e^^'^ — |k|e*'^'''. Here time has been scaled 
by uj so that the dimensionless two-frequency acceleration is 



G{t) = Go - [/m cos(mr) + cos(nr + 0)] (23) 
— Go — f [cos(x) cos (mr) -|- sin(x) cos(nT + 0)] . 

The damping number (7), capillarity number (Fq), gravity number (Go), and 
dimensionless accelerations {fm and /„) are related to the forcing function (23) 
and the fluid parameters by 

_ 2z/A;2 „ _ Ffc3 _ 9ok „ _ gmk . _ gj: 

7= ) i = 5-; tro = — 5-, Jm = :r, In = — ^- l^^J 

Oi) puj"^ uj-^ uj-^ 
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Here v is the kinematic viscosity, F is the surface tension, p is the fluid den- 
sity, and the wave number k is chosen to satisfy the gravity-capillary wave 
dispersion relation 

9ok-r — = (— ) . (25) 



Note that (24) and (25) imply a relationship between the gravity number and 
the capillarity number, namely 



2 



Go + To^— . (26) 



Following [24] , we express the governing equations in the following alternative 
form. We apply {d-,- — 7V^) to (19) to obtain 

{dr - ^V^h - {dr - 7V2) V^={dr- 7V2)J^. (27) 



In (27), we substitute for (dr — 7V^)$ by using (20). The resulting equation 
and equation (19) (which wc rearrange) constitute the system of equations 
that we use in our perturbation analysis: 



FoV' - G{r)] } h = {dr - 7V')J^ + -30 



(28) 
(29) 



3.2 Outline of the calculation 



Wc calculate from the Zhang- Virials equations (28) - (29) systems of ODEs 
for travelling waves in one spatial dimension that are valid for the different 
cases of spatiotemporal resonance described in section 2.2. Our perturbation 
calculations are performed for small amplitude waves in the limit of weak 
damping and forcing (7, /m, /n < 1)- 

For our calculations, we focus on counter-propagating travelling waves having 
critical wave number k (to be determined) which are assumed to be subhar- 
monic to the dominant forcing component cos(mT) and thus, to leading order, 
have frequency m/2. We refer to these waves as the "basic waves". Further- 
more, we insist that /„ not exceed the critical value at which standing waves 
of dominant frequency n/2 bifurcate. Thus, we do not include these waves 
in our nonlinear calculation, and our nonlinear analysis is restricted to the 
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parameter region 

o</^«i, o</„</r* «i 



(30) 



which is bounded away from the bicritical point. 

To facihtate our analysis, we perform four separate calculations, each of which 
pertains to a different possible case of spatiotemporal resonance. The first case 
is a nonresonant case, for which we retain only the basic waves in the leading 
order solution of the perturbation problem. Then we consider cases in which 
the basic waves are (nearly) temporally resonant with their spatial harmon- 
ics having wave number 2k. Based on the arguments in section 2, the three 
resonant frequencies we consider are m (1:2 temporal resonance), |m — n| 
(difference frequency resonance) and m + n (sum frequency resonance). Each 
resonance may be achieved by choosing a particular value of the capillarity 
parameter Fq. For each of these resonant cases, we retain both the basic waves 
and the resonant waves at leading order. For all cases, resonant terms at sub- 
sequent orders in the perturbation calculation lead to solvability conditions, 
from which we obtain travelling wave amplitude equations, and by a further 
center manifold reduction, standing wave amplitude equations. 

3.3 No resonance 

We use the following scahng: 

7 = 671, fn^efl f^^efl + e^fl + ... (31) 
/c = /co + e^/c2 + • • • , dr^dT-\- edr^ + e^dr^ + e^dr^ + ... 
h^ehi + e^h2 + e^h^ + . . . , $ = e$i + e^$2 + e^^3 + ■ ■ ■ 

where < e ^ 1. The fields h and $ are functions of the spatial variable x, 
the fast time r, and the slow times Tj. 

The expressions for and k indicate expansions of the critical wave number 
and forcing value. We find that terms proportional to in and e in A; are 
not necessary. The wave number, forcing, time derivative, and the two fields 
are expanded through 0{e^) because we carry out the perturbation calculation 
to 0{e^). This higher order calculation is needed since A, the cubic coefficient 
in the standing wave equation (18), turns out to be an 0{e) quantity. (This is 
related to a weakly broken time reversal symmetry as discussed in [26].) 

At 0(e), (28) is 

£o/ii = (32) 
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where 



Co = dl + V{Go-Todl). (33) 



Equation (32) has an infinite-dimensional solution space consisting of all plane 
waves Q^^o^+'^^i.^o)'^ that satisfy the dispersion relation 

Q\ko) = G^ko + Tokl. (34) 



Thus, hi should consist of a superposition of these plane waves, i.e. hi = 
J2ko ^i^o) J where the wave number ko may be any wave number 

that fits into our periodic domain. However, at C(e^), all of the amplitudes 
z{ko) are damped on the slow time scales, except for the case Jl(/co) — "^/2, 
ko — 1. Using this a posteriori justification, we choose hi to include only those 
solutions which may grow on the slow time scales. Therefore, hi consists of 
one set of counter-propagating waves: 

hi = zi e^'^^+^f ^ +wi e^'^^-^f ^ + c.c. (35) 



where Zi and Wi are functions of 7^, 7^, and %. 

At C(e^), C(e^) and 0{e'^) we apply solvability conditions which yield the 
respective equations 



1^ = -71^1 + iViWi (36) 
dTi 

^ ■iiy2Zi + ici\zifzi + ic2\wi\'^Zi (37) 



^ = -73-^1 + ^^ys'W^i + csl^ip^i + Cilwil'^zi (38) 



+ ic^lwifzi + ice\wifwi + icjzfwi 



and similar equations for wi which are related by the spatial reflection sym- 
metry 

X ^ -X : {zi,wi) ^ {wi,zi). (39) 



The coefficients in (36) - (38) are given in the appendix. 

We reconstitute the time derivative and amplitudes in the travelling wave 
equations by multiplying (36), (37), and (38) by e^, e^, and respectively, 
adding the results, and letting ezi — >• zi, ewi — > Wi, and edr^ + e^f?r2 + £^^7-3 —>■ 
dr- We obtain 
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dzi 

(-7 + iiy)zi + ir]Wi + {di + ici)\zi\'^zi + (^2 + ic2)\wi\'^ Zi (40) 
+ id^lzil'^Wi + i(i4|-u;ipu'i + id^zlwi 



and a similar equation for related by (39). The coefficients in the recon- 
stituted equations are 

7 = e7i + e^73, = e^i^2, ?7 = e?7i + e^?73, c^i = ecg (41) 
(^2 — ec4, (^3 = ec5, (^4 = ece, ^5 = ecy. 

(Note that €71 = 7, etc. so that the e's drop out of the final equation.) 

Now we solve for k2, the correction to the critical wave number, and and 
in (31), the forcing amplitudes associated with onset. The condition for 
neutral stabihty of the fiat state follows from (40) and is 

rj^^f + iy\ (42) 



At leading order, O(e^), we solve (42) for to find that 

fin = 2m7i. (43) 

At C(e^), we find that 

^,_ kl{m^ + 8Tor ^ k2{mm' + 8To) 
^ 16m7i 4m7i(n^ — m^) 

7i(/^)^(7nW + n^-4m^) ^ 71/^2 (7m^ - SFp) 
2m'n?{'n? — m?Y 4m 

4m7i(n^— m^)^ 4m 
which is minimized at 

_ 27K7m^-8ro) 2(/i)^ 
' (m2 + 8ro)2 (m2 + 8ro)(n2-m2) ^ ^ 



and so 



3 ^ 7?(29m^ + Wro + 320rg) 
■'"^ 2m(m2 + 8ro)2 ^ ' 

2-fi{f^fm{m* + Sm^Fo - Wn^To - 2n^) 
'n?{m? + 8ro)(n^ — rn?y 
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To reconstitute the expressions for critical forcing and wave number, we recall 
(31) to obtain 



/^^2^.- ^'p^;+'^":^^;;f«rg' (47, 

Zm{m'^ + 81 o) 
^ 27(/„)2m(m^ + Sm^Fo - IGn^Fo - 2n*) 



TT?{m? + 8ro)(n^ — rn?Y 
27^(7m^-8ro) 2{Uf 
(m2 + 8ro)2 (m2 + 8ro)(n2-m2)' ^ ^ 

The superscript c indicates that fm has been set to its critical value. 

Due to the restrictions (30) which we placed on the forcing amplitudes, the 
calculation performed above gives us information only about one side of the 
linear stability boundary. In order to obtain expressions for the entire linear 
stability boundary in - space, we perform a similar linear calculation for 
the case that the dominant forcing component is cos(nr). The critical forcing 
and wave number in this case are 



Fn-'^inK (49) 
_ l{fmfnkl{-?,n'^ + %v?GQkn + 4m2n2 - l^k^m'^GQ + 2m^) 
2m2(G'o + 3roA;2)(n2-m2)2 

7=^(-n2 + 4G'oA;„)(53n4 - ITGn^Gofcn + 'ilQGlkl) 
32ron(12G'oroA;2 - AGl - <dV^kW) 
K - I Umfkl ^'kl{Gok^ + 3klro-2n') 

" ^ 2(Go + 3A;2ro)(n2 - m^) ^ 2(6Goro^^ + ^2 + 9r2^4) • ^ ) 

Here kn satisfies 

(51) 

where ft{k) represents the natural frequency given by the dispersion relation 

n\k) = Gok + rok^ (52) 



These linear results are discussed in section 4.1. 

We now return to the case that the cosimr) forcing dominates and continue 
our calculation in order to determine the cubic coefficient. The critical eigen- 
vector V = (fi,f2) for the travelling wave system (40) is defined up to a 
complex constant, which we choose such that Vi and V2 are complex conjugate 
quantities with real part equal to one. (We choose this particular eigenvector 
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to be consistent with the numerical computation of A, which was performed 
in [25].) We reduce (40) to the critical mode to obtain the standing wave 
amplitude equation 

— + ^nonres\Zl\ ^1- 

(53) 



The cubic coefficient Anonres is calculated through its leading order, namely 
C(e). We find 

■^nonres = 2e [ C3 + C4 — C5 — Cq + C7 ^ ^ ] . (54) 

V 71 / 



We substitute for 1/2 and ci . . . cy in (54) and reconstitute to obtain 



_ -37(5m2 + 2ro) 1817771^ _ 2877^^ 

Anonres - + 10(^2 + 8ro) " + UFq ^ ' 



^ 377771 167777 



5(7772 _ i2ro) (7772 - i2ro)2 



Note that Anonres diverges as Fq w? /12. This divergence reflects the fact 
that the second spatial harmonic of the critical mode is resonantly excited 
when Fq = jVl. We perform the necessary calculation for this case next. 



S.J^ 1:2 spatiotemporal resonance 



Now we perform a calculation to handle the case of resonance involving the 
temporal harmonic; this resonance occurs when the spatial harmonic of the 
basic waves oscillates with frequency 777. The condition is 

n'^{2ko)=m^ (56) 



where Q is given by the dispersion relation (52). Solving (56) for Fq, we see 
that the 1:2 spatiotemporal resonance occurs for Fq = Fi:2 = m'^/12, which is 
the value of Fq at which the nonresonant calculation in section 3.3 diverges. 

The analysis here is similar to that of section 3.3, except that we now include 
the resonant mode in our calculation. Thus, 

hi = zi e^'^^+^f ^ +wi 6^^=^-^^^ +Z3 6^^'=^+*"*^ +W3 e'^^'^-i^-' + C.C. (57) 
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Additionally, since we are interested only in the parameter region near the 
resonance, we expand around the resonant value of the capillarity number: 

To = ri:2 + ef (58) 



A solvabihty condition at 0{e^) yields 



^ = -71^1 + iViWi + ieiZiZ3 (59) 
oil 

= (-74 + ^2^4)2:3 + ie2Zi 

and equations related by a spatial reflection symmetry similar to (39). The 
coefficients are given in the appendix. 

The leading order term in the standing wave cubic coefficient Ai.2 depends 
only on terms in (59) and thus may be determined without carrying the per- 
turbation calculation any further. A reduction of (59) to the standing wave 
equation (18) reveals that the leading term in Ai.2 is an C(e~^) quantity given 

by 

2616274 

Al:2 = 2N- (60) 

e(74 + i^4 



We substitute for the coefficients to obtain 



9(ro-ri:2)2 + i67W ^ ' 



which is vahd for Fq sufficiently close to ri:2. 



3.5 Difference frequency resonance 



Now we perform a calculation to handle the case of resonance involving the 
difference frequency mode, which occurs when the spatial harmonic of the 
basic waves oscillates with frequency \m — n\. This condition may be written 

as 

n'^{2ko) = {m-nf (62) 
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where Q is given by the dispersion relation (52). Solving (62) for Fq, we see 
that the difference frequency resonance occurs for 



To = ^diff = - ^nm + j^rn^- (63) 



The calculation is similar to that of the previous section. We let 

To = ^diff + efdiff- (64) 

Now, hi is given by 

hi = zi e^*^^+^f ^ +wi e^*^^-^^^ (65) 

The solvability conditions at O(e^) and 0{e^) yield 

dzi 



dTi 
dz3 
dTi 
dzi 

Wo 



-7i + zi + irjiWi (66) 
(-74 + ii^i)zz 



2 



■ iv2Zi + iriziz^ + ici\zi\^Zi + ic2\wi\^Zi 

+ ic8|-23p-2i +ic9|«;3p2;i 



dz: 



d% 



^ ^ii^bZs + 'ir2zf + icio\zi\'^zs + icu\wifz3 



+ iculzsl^zs + icnlwsl'^zs. 



and equations for wi and Ws related by a spatial reflection symmetry. The 
coefficients are given in the appendix. The values for 1/2, Ci, and C2 are given 
by (A.2), (A.3), and (A.4) evaluated at Tq = ^diff- 

It is not necessary to carry the perturbation calculation further to determine 
the standing wave cubic coefficient A^iff at leading order. The coefficient A^iff 
has two types of contributions. One type is unrelated to the slaved difference 
frequency mode and is equal to Anonres evaluated at Tq = T^iff. The other 
type is due to the quadratic terms in (66) and results from the slaving of the 
damped difference frequency mode. We find that 

^diff -^nonres (ro = r,,;;) + A,,ff (67) 
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where 



2 2e74rir2 



We substitute for the coefficients to obtain 

n2 [1672(71 - m)2 + 9(ro - Vdifff] {n -m){n- 2mf ' ^ ^ 



3. 6 Sum frequency resonance 

The condition for the sum frequency resonance is 

Q2(2A;o) = (m + n)^ (70) 

where Q is given by the dispersion relation (52). Solving (70) for Fq, we see 
that the sum frequency resonance occurs for 

To = ^sum = + \nm + (71) 

The calculation is almost identical to that of the previous section, and the 
result may be obtained by letting m —m in (69). We find 

■^sum ^nonres (Tq ^ sum) ~l~ ^sum (72) 



where 



n2 [1672(71 + m)2 + 9(ro- r,„^)2](n + m)(n + 2m)2- ^ > 



4 Results 



4-1 Linear results 



We now discuss results that apply to the linear instability of the trivial solution 
of the Faraday wave problem, i.e. the flat interface state. Figures 2 and 3 
contain sample results for the case 7n = 4, = 9, Fq = 2, and various values 
of the damping parameter 7. The data are computed both numerically, using 



20 



8 1 1 1 1 

7=0.2 



6 - 

f 

n 

^ 7=0.1 



7=0.05 



7=0.01 
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Fig. 2. Linear stability boundary in /„i - fn space, the parameter space of the two 
acceleration amplitudes in (23). For a given value of the damping parameter 7, the 
flat interface state is unstable above and to the right of the corresponding curve. 
Dotted lines are numerical data; solid lines correspond to the analytical expressions 
(47) and (49). The two are distinguishable on this graph only for 7 = 0.2. The other 
parameters are m = 4, = 9, ^ = in (23), and Fq = 2 in (28) - (29). 

the method in [25], and from the analytical expressions in (47) - (48) and 
(49) - (50). Figure 2 shows the linear stability boundary in - /„ space. 
Figure 3 shows the critical wave number as a function of the quantity %. Note 
that increasing x corresponds to marching counterclockwise around the linear 
stability boundary of figure 2. 

The expressions for critical acceleration and wave number were derived in sec- 
tion 3.3 by performing a perturbation expansion on the Zhang- Viiials equa- 
tions (28) - (29) for small amplitude waves and weak damping and forcing. 
For arbitrary damping and forcing, the linearization of (28) - (29) is a damped 
Mathieu equation for each Fourier mode pfe(T)e''^^: 

pl + 2^k% + {j'k' + n\k)}p, (74) 
^ k{fmCos{mT) + /nCos(nr + (p)}pk- 
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0.98 



0.96 



0.94 




Fig. 3. Critical wave number kc as a function of Xj shown here for x < Xbc- The 
critical wave number decreases as the bicritical point is approached. Dotted lines 
are numerical data; solid lines correspond to the analytical expression (48). The 
parameters used are the same as those in Figure 2. 

where the natural frequency fl{k) satisfies the dispersion relation (52). Thus, 
we may compare the results (47) - (48) and (49) - (50) with known results for 
the Mathieu equation with weak damping and forcing; see [28], for example. 
Here we focus on (47) and (48) which apply when the bifurcation is due to the 
cos(mr) forcing. This bifurcation corresponds to crossing through the right 
side of the linear stability region. Similar statements hold for crossing through 
the top of the linear stability region, when the bifurcation is due to the cos(nT) 
forcing, in which case (49) - (50) are the relevant quantities. 



At leading order, the critical forcing (47) is proportional to the damping 7. 
There are two correction terms. One correction term is proportional to 7^ and 
is independent of /„. This term always has an overall negative sign and hence 
lowers The other correction term is proportional to 7(/n)^ and is due to 
the second forcing component. The overall sign of this term is determined by 
the quantity 



(75) 
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If s < 0, the j{fnY term has an overah negative sign, and thus the second 
forcing component cos(nr) is destabihzing; that is to say, it pushes the bifur- 
cation to occur at a smaller value of However, if s > 0, the second forcing 
component actually stabilizes the flat fluid surface beyond those values of fm 
where it would have otherwise gone unstable. 

By analyzing the expression for s, remembering that Fq is restricted to the 
interval < Fq < we see that there are three possible cases: 

(1) If m/n < v^, the second frequency component is destabihzing for all 
values of Fq. 

(2) If < m/n < y/l + |\/I0, the second frequency component is stabiliz- 



ing for Fo < F, - "^^ 



(3) If m/n > y I + |\/I0, the second frequency component is stabilizing for 
all values of Fq. 

In short, the secondary forcing component is stabilizing if it is at sufficiently 
low frequency compared to the dominant forcing component. (However, since 
our results apply to weak damping and forcing, the effect of the 7/^ term is 
quite small.) 

The bicritical point Xbc may be determined from (47) and (49), the expressions 
for critical forcing. To leading order, it is given by the simple expression 

Xbc = arctan f j (76) 



where kn is determined by the dispersion relation (51). Note that to leading 
order, Xbc depends on m, n, and the capillarity number Fq, and is independent 
of damping and forcing. Using the bounds on kn that are set by the dispersion 
relation, we see that for a given ratio m/n, Xbc takes on extreme values of 



xlc — arctan ( — j at Fq = (gravity waves) (77) 

v m J 

(n \ ^1"^ 
— j at Fq = m^/4 (capillary waves). 

For m < n, xlc ^ maximum and xlc ^ minimum; the reverse is true for 
m > n. As Fq is changed, Xbc varies smoothly and monotonically between the 
two extrema. Examples are shown in figure 4 for m = 4 and various values of 
n. 

The critical wave number, to leading order, is one. This is simply the dimen- 
sionless wave number determined by the dispersion relation fl'^{k) = (m/2)^. 
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Fig. 4. Bicritical point Xbc (in degrees) versus capillarity parameter Fq for m = 4. 
Lines correspond to the expression in (76). Symbols correspond to a numerical 
computation with 7 = 0.1 ("o") and 7 = 0.4 ("x"). 



where VL{k) is given by (52). There are two correction terms. One is propor- 
tional to 7^ and has an overall negative sign. The other is proportional to 
The overall sign of this term is given by the sign oi m — n. Therefore, the 
presence of the second forcing component shifts the wave number in such a 
way as to "repel" it from the other instability associated with the bicritical 
point. This effect can be seen in figure 3, which shows the critical wave number 
for X < Xbc- 

Finally, we discuss the fast-time dependence of the Faraday-unstable mode, 
which we write as p(t). As demonstrated in [5], p{t) will be harmonic or 
subharmonic to the period of the forcing function (23), namely 27r. Previous 
work has depended on a numerically determined (truncated) Fourier series at 
some point in the linear or nonlinear analysis. For instance, in [5,24,25,29] the 
time dependence of the critical mode is written as 



N 

E 

j=-N 



(78) 
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J 

Fig. 5. Magnitudes of the nine most significant fast-time frequency components in 
a neutrally stable Faraday mode near the bicritical point. The vertical axis (note 
log scale) shows the magnitude \aj\ of the frequency component c*-'^, normalized 
so that the largest component has magnitude one. The horizontal axis shows the 
Fourier index j. The components have been arranged in decreasing order of their 
magnitude. Squares correspond to data from a numerical computation. Circles follow 
from the perturbation analysis in section 3. The parameters used are m = 4, n = 9, 
/„ = 3.61, and = in (23), and 7 = 0.1 and Fq = 2 in (28) - (29). These data are 
for the critical mode, with wave number k; the spatially resonant mode with wave 
number 2k will be dominated by a different frequency component determined by 
the dispersion relation (52) (e.g. |m — n| if Fq = Taiff). 

for the harmonic case, or 



iV-l 
j=-N 



(79) 



for the subharmonic case. Then, the coefficients aj are determined numerically. 
This method assumes no a priori information about the relative importance 
of the frequency components kept in the expansion. 

Our analysis determines the relative importance of the frequency components 
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in p(t) for arbitrary m and n in (23). For our perturbation expansion in 
section 3 we assumed that at leading order, the Faraday waves have frequency 
|m. At second order in the expansion we captured the frequency components 
\n — |m|, n + |m, and |m. At third order, we captured the components 
\n — |m|, \2n — |m|, n + |m, 2n + |m, and |m. 

These resuhs arc consistent with what we find numerically. Figure 5 shows the 
nine temporal Fourier coefficients \aj\ that are largest, arranged in decreas- 
ing order of their magnitude. Note that even near the bicritical point, where 
this data was obtained, the |m frequency component is at least an order of 
magnitude larger than any of the other components. 

4-2 Nonlinear results 

4-2.1 One spatial dimension 

In this section we discuss the nonlinear results of section 3 for the cubic coeffi- 
cient A in (18). We have checked our perturbation results with numerical com- 
putations. Figure 6 shows a sample result of A versus the capillarity parameter 
Fo for m/n = 4/9. The solid line corresponds to an expression which matches 
Anonres (55) to A^iff in (67) and thus is valid for all values of Fq away from 
the 1:2 resonance. This expression diverges at Fq = F1.2 = m^/12 as discussed 
in section 3. The dotted line corresponds to the expression Ai.2 in (61). Addi- 
tionally, we have calculated the relative error in the perturbation results for A 
as a function of the damping 7. For instance, for m/n = 4/9 and % = 75° < Xbc, 
as 7 is varied from 0.05 to 0.25, the relative error in A(Fo = ^diff) increases 
from 0.001 to 0.25, and the relative error in ^(Fq = Fi:2) increases from 0.05 
to 0.43. 

■^nonres, the value of the cubic coefficient away from the 1:2, difference, and 
sum frequency resonances, was computed in section 3.3 and is given by (55). 
Anonres IS proportioual to the damping parameter 7. Furthermore, Anonres is 
always negative indicating that in the nonresonant regime, the bifurcation 
from the fiat state is always supercritical. 

Ai:2, the value of the cubic coefficient near the 1:2 temporal resonance, is 
given by (61). This quantity was derived for |Fo — '^di^f \ — (^(7)- Thus, in 
the region of validity, A1.2 — 0(7^^). This is significantly larger in magnitude 
than A^onresi which is ^(7). Furthermore, A^.^ is negative, again indicating a 
supercritical bifurcation. This large negative contribution is manifest as the 
large dip around Fq ~ Fi:2 = 4/3 in figure 6. A\.2 has a global minimum 

2 

at (Fq, Ai:2) = (Fi:2,— ^), SO that exactly at the 1:2 resonance, the value 
of the coefficient A is inversely proportional to the damping, as predicted by 
the symmetry arguments of section 2.2. Thus, near the 1:2 resonance, one- 
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Fig. 6. Cubic coefficient A in (18) as a function of the capillarity parameter Fq. 

The dots correspond to a numerical computation. The dotted line corresponds to 
the expression for Ai.2 in (61). The solid line corresponds to an expression which 
asymptotically matches Anonres and A^iff (details not given). The large dip at 
Fq ~ Fi:2 = 4/3 is due to the 1:2 resonance discussed in section 2.2. The small 
bump around Fq ~ ^diff = 17/6, at which the one-dimensional waves have their 
largest amplitude, is due to the difference frequency resonance, also discussed in 
section 2.2. The parameters used are m = 4, n = 9, % = 75° and ^ = in (23), and 
7 = 0.05 in (28) - (29). 

dimensional waves will decrease significantly in amplitude. 

Adiff, the value of the cubic coefficient near the difference frequency resonance, 
is given by (67). The condition for difference frequency resonance is Fq = ^diff, 
where T'diff is given by (63). Since Fq is restricted to the range [0,m^/4], this 
condition can only be met for certain m/n ratios. Specifically, r^iff G [0, m^/4] 
only for 

m/n e Ml U Ma, Mi = [^2 - 1, 2 - ^2], M2 = [2 + v^, oo). (80) 

Thus, while the 1:2 resonance was relevant for all possible forcing frequency 
ratios m/n, this is not the case for the difference frequency resonance. The 
difference frequency resonance results in a contribution to A, namely A^iff 
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given by (69), and thus A has a local extremum at Fq = ^diff- The sign of 
Adiff is given by the sign of n — m. If the secondary forcing component is at 
a higher frequency than the primary, i.e. if m/n e Mi, then the difference 
frequency resonance results in a positive contribution to A. The extremum is 
a local maximum, and the amplitude of the supercritical waves increases as 
the resonance is approached. This is demonstrated by the small bump around 
To ~ Tdi// = 17/6 in figure 6. If m/n G M2, then the contribution is negative. 
The extremum is a local minimum, and the amplitude of the waves decreases. 
In either case, the extra contribution to AITq = Taiff) that is due to the 
difference frequency resonance is proportional to {fnY/l ^ predicted by the 
symmetry arguments of section 2.2, and thus is a smaller affect than the 1:2 
resonance. 

For the case that m/n e Mi, when A has a local maximum, it is possible 
for this maximum to actually cross the A — Q axis and become positive, 
thus causing the bifurcation to become subcritical. This will only happen if 
/„ is sufficiently large. Writing A^iff as adif}{fnf'-i then the condition for a 
subcritical bifurcation is 



fn > 



Of course, condition (81) must be considered in conjunction with the condition 

fn < frv which we enforced in (30). (An example of a subcritical bifurcation 
may be obtained with the parameters m/n = 49/100, Fq = ^diff = 2801/12 fa 
233.4, 7i = 0.01 and = 3.9 < 3.94 in which case A = 0.57 > 0.) 

Now we turn to the results for the sum frequency resonance. Asum is given by 
(72). The condition for sum frequency resonance is Fq = Tsum, where Tsum is 
given by (71). Similar to the difference frequency resonance case, this condition 
will only be met for certain m/n ratios. Specifically, the sum frequency mode 
resonance is possible only for 

m/n >V2 + 1. (82) 



Thus, the sum frequency resonance can only be realized when the second 
forcing component is at sufficiently low frequency. The sum frequency reso- 
nance results in a contribution to A, namely Agum, which is given by (73). 
Agum has a local extremum at Fq = F^^^. Like the difference frequency case, 
the contribution to A due to the sum frequency resonance is proportional to 
ifuY/l- Unlike the difference frequency case, the Agu,,, contribution always 
has a positive sign. However, this contribution is generally not significant be- 
cause the algebraic pref actor in A^^^ is small for values of m/n for which the 
sum frequency resonance is possible. 
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1)2 C-) 1:2 (-) hli-l mi-} 1:2 (-) 
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^, I I ^ smiH I mmM ^ 

m/n = V2-I 2-V2 V2 + 1 2 + V2 

Fig. 7. Regions of forcing frequency ratio m/n in which the 1:2, difference, and 
sum frequency resonances are possible for one-dimensional waves. The plus (+) and 
minus (— ) signs indicate whether the resonance results in a positive or negative 
contribution to the cubic coefficient A in (18). In the case of (— ) the bifurcation 
to one-dimensional waves is necessarily supercritical. Note that only points corre- 
sponding to rational numbers on the number line are meaningful. 

A partial summary of the results for 1-d resonances may be found in figure 7. 
This number line shows the regions of forcing frequency ratio m/n in which 
each type of resonance is possible. The plus {+) and minus (— ) signs indicate 
whether the resonance results in a positive or negative contribution to the 
coefficient A, and hence whether it makes the supercritical waves larger (-I-) 
or smaller (— ) in amplitude. 



4-2.2 Two spatial dimensions 

We now present nonlinear results for Faraday waves in tw'o spatial dimensions. 
We have computed the cross-coupling coefficient B{6) in (8) using the method 
in [25]. We interpret features of B{6) in light of the resonances discussed in 
section 2.2. Many of these features may be understood by means of a simple 
argument which is valid for weak damping and forcing. We simply solve the 
spatial resonance condition (2) or (4) for 9i:2, Osum or Odiff, which are the 
angles at which the 1:2, sum frequency, and difference frequency resonances 
occur. To do this, we must set = k{Q) where k{Q) is the inverse of the 
dispersion relation (52) and fl = m, m + n oi \m — n\ depending on the 
resonance under consideration. A number of results immediately follow: 

• The 1:2 resonance is possible only for Fq > m^/12 = F1.2. 

• The difference frequency resonance is possible only for 

m - \\l2m? + 24Fo < n < m + \\l2m? + 24Fo. (83) 

• The sum frequency resonance is possible only for n < —m + + 24Fo. 
Prom these statements, we also see that 

• The ranges of 6*1:2, Qsum and Qdijf are restricted. 

• There are some forcing frequency ratios mjn for which the sum and differ- 
ence frequency resonances are not possible for any value Fq. 
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Fig. 8. Cross-coupling coefficient B{9) in (8). The solid line corresponds to 
m/n = 8/9 in (23); the dotted and dashed lines correspond to m/n = 8/11 and 

8/21 respectively. For each curve, the parameter x is chosen to obtain a harmonic 
instability near the bicritical point. The other parameters are ^ = in (23), and 
To = 14 and 7 = 0.1 in (28) - (29). The large dip at 6* = 01:2 ~ 70° is due to the 1:2 
temporal resonance discussed in section 2.2 and is independent of the second forcing 
component. The small spike is due to the difference frequency resonance discussed 
in section 2.2. We have removed from this plot the region near 6 = 60° where B{9) 
necessarily diverges; a calculation for the hexagonal lattice is required here. 

An example is given in figure 8, which shows the cross-coupling coefficient 
B(9) computed for forcing frequency ratios m/n = 8/9, 8/11 and 8/21 for 
fixed fluid parameters; % is chosen in each case to obtain a harmonic instability 
near the bicritical point. The large dip at ^ = 9i:2 ~ 70° is consequence of 
the 1:2 resonance. At this angle, which is predicted by the weak damping 
argument given above, there is a resonant triad comprised of two Faraday- 
unstable modes with dominant frequency m/2 and the weakly damped mode 
oscillating primarily with the harmonic frequency m. As expected from the 
analysis of section 2.2, near this angle, the weakly damped mode contributes 
to B{9), which here is manifest as the large dip. This phenomenon is similar to 
the 1:2 resonance in one spatial dimension, which resulted in a large dip in the 
cubic self-interaction coefficient A in (18). It follows from the weak damping 
argument that 6*1:2 will depend on m and Fq but will be largely independent 
of the parameters n, /„, and 7. The independence with respect to n is evident 
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in figure 8, in which the dip occurs at the same angle for m/n — 8/9, 8/11, 
and 8/21. 

For all numerical calculations that we performed, the 1:2 resonance resulted 
in a large negative contribution to the cross-coupling coefficient B{9). As dis- 
cussed in section 2.1. this type of contribution is destabilizing for superlatticc 
patterns with characteristic angles 6'/inear 9i.2. Our numerical results (not 
shown) indicate that the magnitude of the dip caused by the 1:2 resonance 
follows the scahng law that we deduced from symmetry considerations in sec- 
tion 2.2, and that we derived for one-dimensional waves: namely, that the 
contribution from the weakly damped mode scales like I/7. 

The sum frequency resonance angle 9sum ^lay also be predicted by the weak 
damping argument. However, unlike the 1:2 resonance described above and the 
difference frequency resonance described below, the sum frequency resonance 
for two dimensional waves is quite difficult to detect numerically for typical 
values of m/n and for small 7. This is consistent with the result for one spatial 
dimension, and consistent with the fact that the mode oscillating at the sum 
frequency has a larger wave number and thus is more strongly damped. 

Finally, we turn to results for the difference frequency resonance. The effect 
of the difference frequency resonance may be seen in figure 8, and is manifest 
as a spike in the plot of B{9). Let us first concentrate on the solid curve in 
figure 8, which corresponds to a forcing frequency ratio oi m/n = 8/9. For 
this case, at 9 — 9diff 17° there is a resonant triad composed of two modes 
with dominant frequency m/2 and the weakly damped mode oscillating with 
dominant frequency \n — m\. As expected from the analysis of section 2.2, near 
this angle, the weakly damped mode contributes to B{9), which causes the 
spike. This phenomenon is similar to the difference frequency resonance in one 
spatial dimension, which resulted in a contribution to the cubic self-interaction 
coefficient A. 

As with the case of 1:2 resonance, the simple argument we have used to predict 
the resonance angle 9diff relies only on the dispersion relation (52) and on the 
trigonometric relations (3) and (5). By examining these two expressions, we 
expect that 9diff will depend on m, n, and Fq but will be largely independent 
of the parameters /„ and 7. The dependence on the second forcing frequency 
n is evident in figure 8, in which shifting from n = 9ton = ll causes the 
spike to shift from 9diff ~ 17° to 9diff ~ 47°. 

Figure 9 shows the angle of spatial resonance 9diff versus the capillarity pa- 
rameter Fo for the forcing frequency ratios m/n = 8/9 and 8/11 and for 
various values of 7. The solid lines represent the prediction of 9diff based 
on the weak-damping argument, while the points represent data from a full 
numerical computation of B{9). 
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Fig. 9. Angle of spatial resonance Odiff versus capillarity number Tq. Lines cor- 
respond to a prediction of Odiff based on the dispersion relation (52) and on the 
trigonometric relation (5). Symbols correspond to a numerical calculation of Odiff'- 
7 = 0.2 ("x"), 7 = 0.8 ("□"). The other parameters are m = 8, x = 50° and = 
in (23). 

Another result that follows from the weak damping arguments is that if the 
second forcing frequency n is sufficiently different from m, the difference fre- 
quency resonance will not be possible for any value of Fq. This phenomenon is 
demonstrated in figure 8. The forcing frequency ratio m/n = 8/21 violates the 
condition (83) for all allowed Fq, and the corresponding B{9) curve (dashed 
fine) displays only the 1:2 resonance effect. 

Now we discuss the magnitude and direction of the difference frequency mode 
resonance effect. In contrast to the 1:2 resonance, we find that the difference 
frequency resonance may result in a spike or a dip. Limited numerical results 
for the sign of the resonance effect agree with the result for one spatial dimen- 
sion discussed in section 4.2.1. In particular, we have performed computations 
at 7 = 0.1 for the forcing frequency ratios m/n = 8/7, 8/9, 8/11, 10/7, 10/9 
and 10/11, each for values of Fq ranging between and Tmax = itl'^/4:. In all 
cases, we observe that if n < m then the difference frequency resonance results 
in a dip at 9 — Odiff] H n > m, it results in a spike. An example of the former 
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Fig. 10. Cross-coupling coefficient B{6) for forcing frequency ratio m/n = 8/7 
in (23). The difference frequency resonance results in a small dip around 

6 = ^diff ~ 17°. This is in contrast to the case m/n = 8/9, which produces a 
spike at the same angle. The other parameters are % 37° and (f) = in (23), and 

7 = 0.1 and To = 14 in (28) - (29). 

case is shown in figure 10, which corresponds to a forcing frequency ratio of 
m/n — 8/7. The resonant wave number and the angle of spatial resonance 
Odiff are (approximately) the same as for the case m/n — 8/9 shown in fig- 
ure 8; however, the difference frequency mode resonance now results in a very 
small dip rather than a spike, as before. 

As in the one-dimensional case, the magnitude of the difference frequency res- 
onance effect follows the scaling law that we deduced from symmetry consid- 
erations in section 2.2, namely that the contribution from the weakly damped 
mode scales like {fnY /l- This scaling may be seen in figure 11. We compute the 
magnitude of the effect by finding B{6diff ) - Bf^=o{ddiff), where Bf^=o{Odiff) 
is the value of the cross-coupling coefficient at the resonant angle computed 
without the second forcing component. We plot the size of the spike versus 

(/n)V7- 

As discussed in section 2.1, a spike occurring at spatial angle 9 = Odiff will help 
stabilize SL-I patterns with characteristic angles 6h near Odiff - To demonstrate 
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Fig. 11. B{9^iff) — Bf^^Q^Odiff), the magnitude of the difference frequency spike, 
versus {fnf/l- The damping parameter 7 is varied between 0.01 and 0.1, and the 
strength of the second forcing frequency is varied between (which corresponds to 
single frequency forcing) and Best-fit lines are also shown. The other parameters 
are m = 8 and ^ = in (23), and Tq = 14 in (28) - (29). 

this effect, we consider an example form/n = 8/ll forcing, with 7 = 0.2, and 
To = 13 and focus on the case of a harmonic instabihty. These dimensionless 
parameters can be realized, for instance, by a fluid with surface tension T = 
4.2 dyn/cm, density p = 1.0 g/cm"^, and kinematic viscosity u = 0.01 cm^/s 
being forced with base frequency c<j/(27r) = 16.2 Hz. (The fluid properties here 
are similar to those of water, but with lower surface tension. This situation 
might be achieved by the use of a surfactant). 

When X = 60.5°, there is a spike in B{6) at Odiff = 46.9°, which is close 
to the value of 47.0° that is predicted by the weak-damping argument. We 
have performed a limited bifurcation analysis similar to that in [25]. The 
stability of SL-I patterns is computed with respect to simple rolls, simple 
hexagons, and various rhombic patterns. An SL-I pattern with lattice angle 
6h ~ 47° is stable for a small range of / above, but bounded away from, 
onset. A higher order calculation is necessary to determine whether it is the 
superhexagon or supertriangle variety of SL-I pattern that is stabilized (these 
two different types of SL-I patterns have different phases associated with the 
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Fig. 12. Superhexagon (left) and supertrianglc (right) patterns with characteristic 
angle 9^ ~ 47°. For m/n = 8/11 forcing with 7 = 0.2 and Fq = 13, both patterns are 
unstable for x = 0°. For x = 60.5° < Xhc-, oiie of these patterns is stabilized by the 
difference frequency resonance effect; a higher order calculation is needed to deter- 
mine which one. The patterns shown were created by an appropriate superposition 
of the twelve critical Fourier modes [18]. 

complex amplitude). When x = 0° (i.e. when there is only single frequency 
forcing) the spike in B{9) due to the difference frequency disappears and the 
SL-I patterns with Oh ~ 47° are unstable. Superhexagon and supertriangle 
solutions with Oh ~ 47° are shown in figure 12. 



5 Conclusions 

In this paper we have examined the role that weakly damped modes play in the 
pattern selection process for Faraday waves forced with frequency components 
muj and rvjj. Our symmetry arguments predict that the modes oscillating 
primarily with the harmonic frequency mo;, the difference frequency \n — m\uj, 
and the sum frequency (n + m)a; will be the most important in terms of 
their contribution to the cubic coefficients A and B[0) in the standing wave 
equations (8). The symmetry considerations also provided scaling laws for the 
magnitude of these resonance effects. 

Starting with the Zhang- Viiials Faraday wave equations, we performed a 
weakly nonlinear analysis for weak damping and forcing in order to calculate 
expressions for the self- interaction coefficient A. We obtained expressions for 
the critical forcing acceleration and wave number, and analyzed them to eluci- 
date the role played by the secondary forcing component. We also were able to 
identify the most important frequency components in the unstable eigenmode 
in terms of the integers m and n. We then analyzed the expression for the 
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Fig. 13. Cross-coupling coefficient B{6) in (8) computed for the case of a 
four-frequency forcing function analogous to (23), with frequency components 
(m, n,p, q) = (8, 9, 11, 13). For this calculation, the m = 8 forcing component dom- 
inates, but we are near a "quad-critical" point so that the other three frequency 
components are strong. The three spikes at 6 = 13°, 47°, 73° are due to reso- 
nances with the weakly damped modes oscillating with the difference frequencies 
n — m, p — m and q — m. The other parameters used are Fq = 12.4 and 7 = 0.25 
in (28) - (29). The stability of the SL-I patterns with Oh ~ 13° is enhanced by the 
additional spikes at « 60° ib 9^; see [25] for a detailed discussion of how these two 
additional symmetry-related angles are relevant. 

cubic coefficient A, and determined the sign and scaling of contributions due 
to the various resonance effects. 

Wc then used the Zhang- Vihals equations to numerically calculate the cross- 
coupling coefficient B{6) according to the method in [25]. The predictions of 
the symmetry arguments were manifest. The results for B{6) are of particu- 
lar interest since this coefficient is crucial in determining the stabihty of SL-I 
patterns fike those observed in [2] . We made use of a weak-damping argument 
relying only on the dispersion relation to successfully predict the resonant an- 
gle. While our symmetry arguments predict the scaling of the resonance effects 
and our weak-damping argument predicts the angle, neither argument predicts 
the sign of the contribution to B{9). Our numerical calculation revealed that 
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the 1:2 resonance results in a dip, and thus is destabihzing for SL-I patterns. 
However, the difference frequency resonance in some cases results in a spike, 
which can help stabilize SL-I patterns with characteristic angles near the res- 
onant angle ^di//- This was demonstrated by means of a simple bifurcation 
example. 

We may now speculate on the role of the bicritical point in stabilizing SL-I 
patterns. It is been observed that SL-I patterns occur in experiments only 
for parameters near the bicritical point. It is tempting to believe, then, that 
the weakly damped mode associated with the secondary forcing component is 
somehow responsible for the pattern. Here we have shown that this interpre- 
tation is not necessarily the correct one. Proximity to the bicritical point (i.e. 
making /„ as large as possible before switching over to the other instability) 
maximizes the strength of the difference frequency mode. As we have seen, 
this mode can help stabilize the SL-I pattern. 

It will be interesting to consider Faraday waves forced with more than two fre- 
quency components. In this case, more difference frequency mode resonances 
will be possible. As demonstrated in figure 13, if the parameters are chosen 
properly, multiple spikes in the cross-coupling coefficient B{9) might conspire 
to further enhance the stability of a particular superlattice pattern. 



A Appendix A 



We now give the expressions for the coefficients in the travelling wave equations 
(36) - (38), (59) and (66) which we computed in section 3. 
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